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Recent experimental data on electromagnetic and hadronic associated strangeness production suggest the possible existence 
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We present a non-relativistic quark molecular model (QMM) of the strange crypto-exotic pen- 
taquark baryon spectrum motivated by the recent data showing narrow, resonance enhancements 
in electromagnetic and hadronic production of kaons. Our model assumes color octet bonded quark 
molecular clusters forming exotic color singlet pentaquark baryons. We develop explicit molec- 
ular pentaquark wavefunctions that exhibit color interchange and cluster fermion antisymmetry. 
Our color electro-dynamics (CED) inspired Hamiltonian, which includes confinement, one gluon 
exchange and color magnetic interactions, is a natural generalization of the Isgur and Karl quark 
model Hamiltonian that reproduces the conventional meson (qq) and baryon (qqq) spectrum [1] . We 
introduce the idea of Color Magnetic Confinement (CMC) which is employed to restrict the physical 
pentaquark spectrum and constrain the fundamental constants (color gyromagnetic ratios) of the 
^ ' QMM-CED Hamiltonian. Comparing our QMM spectrum predictions with the recent experimen- 

' tal data in associated strangeness photo- and electro-production reactions, we identify 4 candidate 

. states that have masses and total widths consistent with exotic pentaquark structure. 

^ I PACS number(s); 12.39.Jh; 24.30. Gd; 13.40.-f; 13.75.-n 
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' of narrow multi-quark 5 = [2,3] and S = -1-1 [4,5] states near 2.0 GeV. These pentaquark systems are in the realm of opening 
^— ^ a complete new area of exciting physics on both experimental and theoretical fronts. 

■~>«^ ' The existence and corresponding properties of exotic matter remains a puzzling open question in hadron spectroscopy. 

i^H ' Resolving this mystery is important for understanding how QCD provides a consistent description for the binding of quarks 
I , and gluons into hadrons [6]. Exotic hadrons are color singlet multiquark combinations (qqqq, qqqqq . . . ) beyond the standard 
^ valance qq-(meson) and qqq-(baryon) quark content. Crypto-exotics are those multiquark states which have degenerate quantum 
^ I numbers with the normal hadrons (i.e., non-exotic quantum numbers). Examples include the ao//o(980) resonances, which are 
. established scalar meson states that find their most natural interpretation as hadronic molecules, having wave functions with 
^ ' large components of lightly bound KK mesons [7,8] . The most stable configurations of qqqqq exotics are believed to be strange 
, molecules, like 5* = (uds)-(us) ox S = -1-1 (udd)-(ds), bound by long range color forces acting between color octet quark 
• clusters (octet-bonded states). Exotic molecules such as (qq)-(qqq) can exist in the color 6-6 representations (sextet-bonded 
$_( ' states), however these are expected to be much less stable because a color triplet qqq subsystem always contains a qq pair which 
can form a color singlet, allowing the sextet bonded qqqqq system to easily dissociate into a normal color singlet baryon (qqq) 
and meson (qq) [9,10]. For a baryon B = qqq and a meson M = qq, if B and M are separately color singlets, the molecules 
{B — M) will interact only via the relatively weak Van der Waals forces of QCD and thus will dissociate rapidly (and have 
large decay width) into a normal baryon and meson. However, if the molecular system [B — M) is a color singlet with B and 
M being separately color octet states and possessing non-zero relative orbital angular momentum L, this system should be 
linked by very strong color forces and exhibit enhanced stability (i.e., long lifetime and narrow decay width compared to usual 
hadronic resonance states) since the anti-quark in M is prevented by a centrifugal barrier from tunneling into B where it can 
combine with a quark to form a color singlet meson. Additional stability is anticipated if at least one molecule cluster contains 
a strange (or charm) quark, since then decays to open non-strange {ti,p,uj) channels will be OZI [11] suppressed. 

For S = G strange nucleon resonances, each color octet molecule cluster must have a strange quark, B = qqs and M — qs. 
This naively paradoxical feature of having 5 = with strange exotic structure generates an interesting selection rule for 
detecting exotic baryons: the most stable pentaquark baryons will preferentially decay into strange meson and baryon channels 
by the OZI rule (such as K^A, K*'^A, K^T,, etc.), since the strangeness in each molecular B and M cluster is conserved in the 
OZI enhanced final states. Consequently, the S — pentaquark exotics will potentially show up as narrow nucleon resonances in 
kaon {K or K*) photo- and electro-production reactions, yet totally decouple from the non-strange n,p,Lj production channels 
(or possibly make very small OZI evading/violating contributions). 

This paper is the first work to suggest the possibility of exciting 5 = pentaquark resonances in single kaon electromagnetic 
production from the nucleon (proton or neutron). The recently reported 5 = -1-1 pentaquark baryon resonance in K'^ K~ 
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photo-production from the neutron [4] is a crypto-cxotic state with A quantum numbers (similar to the A (1405) which is widely 
accepted to be a p molecule). In this work we restrict our analysis to the S = nucleon resonance sector, leaving the 
S — ±1 exotics for consideration iu a follow up paper. Before passing into the formalism section, we note that previous authors 
claim that the S = strange exotic baryons and their excited states are expected to fall on linear Regge trajectories [12,13]. 

The excited molecular pentaquarks presumably also have anomalously narrow decay widths in comparison to normal excited 
baryons since the same enhanced stability arguments apply. In contrast, our Quark Molecular Model (QMM) implementation 
permits only 4-low lying S = pentaquark solutions, and dynamically excludes the possibility of additional excited states. 
This limited spectrum feature is a novel prediction of this QMM and stands in contrast to the usual outcome of quark model 
spectroscopy: to generate an overabundance of excited states, and providing theoretical bias for the experimental dilemma of 
"missing resonance^' . 

This paper focuses on two main areas: (1) a dynamical quark model of the 5 = strange exotic baryon (us)-(uds) spectrum; 
and (2) new/recent experimental data that provides preliminary (albeit inconclusive) support for the possible existence of 
strange exotic resonances. 

Section II describes the general formalism of our quark molecular model (QMM), and section III describes the current status 
of experimental data that suggest resonance enhancements due to the S = pentaquark baryons. In forthcoming papers we 
will present a generalized hadrodynamical Lagrangian model that incorporates exotic resonances through strong final state 
interactions in kaon photo- and electro-production reactions [2,3], where we will demonstrate that the exotics have a strong 
forward angle production sensitivity. 



II. FORMALISM 



Our approach is to first construct totally antisymmetric spin-flavor-space-color pentaquark baryon wavef unctions, and then 
calculate the low lying crypto-exotic baryon spectrum using a non relativistic quark molecular model (QMM) basis for the 
spatial wavefunctions and energy eigenvalues. In this work, we assume the pentaquark crypto-exotic systems are heavy, color 
octet bonded {B — M) molecules with anomalous stability. Dynamical stability and explicit decay width calculations will be 
investigated in a future work. Based on the ideas of our introductory discussion, the lightest, stable S — exotic baryons are 
expected to have two distinct features: (1) contain strange quark molecule clusters B = qqs and M = qa to be OZI protected 
from non-strange decay channels [11]; (2) molecule clusters possess at least one unit of relative orbital angular momentum to 
generate a centrifugal barrier which suppresses quaxk rearrangement decay processes. In the derivation that follows, we show 
that both of these qualitative features naturally arise from an explicit minimal construction of pentaquark wavefunctions with 
color octet bonded molecular substructure when total antisymmetry (Fermi statistics) of quark interchange is enforced. 

In order to construct a minimal color singlet pentaquark molecule with total fermion antisymmetry, first consider the normal 
3-q baryon SU{3)f flavor octet and decuplet wavefunctions. Normal ground state 3-q baryons have totally antisymmetric 
quark combinations which are products of color [antisymmetric (A) singlet: Ci], flavor [mixed symmetry (MS) octet: -FmSi or 

symmetric (S) decuplet: Fg ], spin [mixed symmetry (MS) spin doublet: Xms' '^^ symmetric (S) spin | quartet: xj]i and 
spatial [symmetric (S) even parity: wavefunctions, ^ symmetry(color, flavor\spin^'^^'^*^): 

*A(l,8|i^)=Ci (F|fsxL)^*s {octet) (1) 
*a(1, 1o|^) = C\ Fi° x| *J (decouplet), (2) 

where the large round brackets denotes the symmetric product of MS flavor and spin functions. The corresponding anti- 
symmetric (Fx) A product is excluded for ground states {L = 0, and even parity) because the color wavefunction is totally 
antisymmetric under quark interchange, thus the ground state 3-q baryons always have an overall symmetric product of flavor- 
spin-spacc wavefunctions. It is interesting to note that there is no flavor singlet (antisymmetric u,d,s wavefunction) 3-q state 
in the observed low lying baryon spectrum (the A(1405) is not a candidate since it is interpreted as a K^p molecule, and the 
ground state A is the isospin singlet member of the flavor octet). 

In the quark model, it is easy to see why there is no flavor singlet baryon (at least, in the or |^ ground states). 
Constructing the quark model as before, and requiring that all physical hadrons must exist in color singlet configurations, we 
obtain overall MS and S symmetry for the i and | color singlet (g) flavor singlets: 

*M5(l,l|i^) = CiFlxL*J (3) 

^s{l,l\l^) = C\FXxl^t- (4) 

Hence ground state flavor singlet baryons arc inconsistent with the total anti-symmetry requirement. However, if the flavor 
singlet FX{uds) i"*" and |^ baryons have octet color symmetry, then they can bind with qq color octets to form physical 
color singlet baryon states. Our QMM is built on the assumption of color octet bonded a-f3 molecules between baryon number 
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Ba = a(u s) and Bg — 1 I3{uds) basis states. As previously discussed, we restrict our attention in this paper to baryon 
molecules with quantum numbers 5 = 1, isospin 7 = i, and hypercharge Y = B + S = 1, which are crypto-exotics in the 
nucleon resonance spectrum (i.e., same quantum numbers as the N* states). 



A. Color Octet a(u s) and /3(uds) Bsisis States 

We construct totally antisymmetric, color singlet pentaquark baryons which are bound color octet quark molecule clusters 
in the form: (as — /3s)i, whore as and /^g are generic labels referring to the color octet B = a(us,ds) and B = 1 /3(uds) 
clusters respectively. The central idea required to derive explicit wavefunctions in the a-P QMM basis is to exploit the 5*17(3) 
symmetry of usual quark model color and flavor functions to generate pentaquark states that are color singlets and satisfy 
quark interchange antisymmetry. 

Consider a color octet baryon state (which, like colorful quarks, can not exist unbound in isolation) generated by a SU{3) 
Color •<-> Flavor symmetry transformation: 

CA{r,b,g) ^ CMs{r,b,g) = Fms (w -> r, d -» 6, s -> g) (5) 
FMsiu,d,s) ^ FX{u,d,s) = C\{r ^ u,b ^ d,g ^ s), (6) 

where C\ and F*^c,- arc the standard antisymmetric color singlet and mixed symmetry flavor octet functions written in the 
literature describing details of the conventional quark model. The color-flavor (C/F) transformation simply takes the color 
singlet function, replacing labels (r,b,g) with (u,d,s) to generate the flavor singlet flavor function, and the flavor octet functions, 
replacing labels (u,d,s) with (r,b,g) to generate the color octet functions. The C/F transformation generates one color-octet, 
flavor-singlet baryon, ^(i ): 

l^8|(8,l)|i^> =F\(u,d,s) (c^s(r,6,ff)xL)^*s , (7) 
where the MS color and spin function product is symmetrized, and one color-decuplet, flavor-singlet baryon, ^(|^): 

|6o|(10,l)|^^> =F\{u,d,s)Cf{r,h,9)xl^^- (8) 

Our goal is to produce overall antisymmetric pentaquark wavefunctions from color octet bonded molecular clusters a(q s) 
and /3(uds). We accomplish this by forming the anti-symmetrized product of a (8) /? quark cluster wavefunctions with odd 
{L = 1) relative orbital angular momentum: 

*^,3(J+) = iM)^llB = l,/=i F = 1|J+ > (9) 



a{q-s\Q-,l-)\\ Cjtis(a) Fl(a) ^\]s{a) <l>s(a) > (10) 



® I P{uds\^^) II CltsiP) xisW) > 

= C\{aP\X) FiiaP\X) (11) 

Here, we suppress Clebsch-Gordon factors and only express how overall fermion and color antisymmetry is generated in the 
QMM basis. The idea is to simultaneously couple the color octet a and [3 functions to an overall antisymmetric (a — /3)\ color 
singlet function (to satisfy the quark confinement requirement). For notation when we consider experimental data later, we 
refer to the physical [a — (5) pentaquark molecules as the "X" states. In the QMM construction, one can see that these X 
states must have even parity and the B = 1, 1 = = \ quantum numbers of nucleon resonances. The color singlet function 
C\{af3\X) (which couples ?>a ® ^ Ix) is totally antisymmetric under color exchange. We construct explicit antisymmetric 
pentaquark color functions in the next subsection. 

Note the flavor doublet Fg{af3\X) and spin ■)^^^'^^{a(5\X) functions individually have 5-symmetry, and the total X angular 
momentum J = L + Sx couples total spin Sx = Sa + Sp and orbital angular momentum L (with appropriate Clebsch-Gordon 
factors). A detailed analysis shows that the quarks in the Q(qs) cluster can couple to either spin S = 0, 1, in which case 
the coupled spin-space functions for a must be in ^ or S configurations respectively (for Sa ~ 0,1 spin state orthogonality). 
Fermion interchange antisymmetry arises from the product of ^-symmetry color and S-symmetry fiavor, spin, and space {a(3\X) 
functions. The S-symmetry of the iso-doublet molecular flavor function Fg{ap\X) is achieved from the product of two totally 
anti-symmetric flavor functions: the iso-doublet FX{a) and flavor singlet FX{f3). It is important to note here that interchange 
between quarks in |a > and quarks in |/3 > only preserve the overall exact fermion antisymmetry if the (a — /3) molecule has even 
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orbital angular momentum L, because otherwise (for odd L) a ^ [3 quark exchange would generate a MS fermion exchange 
symmetry component due to the negative phase coming from the odd parity orbital (a — 0) spherical harmonic Yj^j^{al3\X) 
function. Our QMM restriction to odd-L states is essentially a dynamical constraint. Primarily, even-L = 0, 2, 4, ... states could 
not be formed in a cluster basis because certain quark exchanges (like s from |/3 > exchanged with w or d from \a >) would 
generate wavefunctions which are not cluster factorized (products of |q; > g) |/3 > cluster functions) and hence destroys the 
chance for defining "basis states" with invariant multiplet quantum numbers. Furthermore, a ^ quark exchange would allow 
s or s tunneling and subsequent annihilation, which obviously destabilizes the pentaquark bound state. In QMM we prohibit 
a <-> /3 quark exchange by requiring the ground state to have one unit of relative orbital angular momentum, and excited states 
must take odd-L (to generate a centrifugal barrier to dynamically suppress quaxk rearrangement decay and to provide the even 
parity required for crypto-exotic nucleon quantum numbers). So the ground state L = 1, and L = odd excited state QMM 
wavefunctions have a restricted cluster fermion antisymmetry (that is, exact A-cxchange symmetry for permutations within the 
a{u,d,s) cluster and/or permutations within the /3{q,s) cluster, but MS-symmetry for interchanges between clusters). Exact 
overall fermion antisymmetry could be generated for a <-» /3 quark exchange in odd-L molecules if mixed color symmetry were 
allowed, but color confinement requires all physical states to be A-co\ot singlets. Prom this more fundamental point of view, 
our QMM cluster fermion antisymmetry is the dynamical symmetry enforcing exact color confinement in a basis of factorized 
quark cluster functions with odd relative orbital angular momentum (and nucleon. A*'* quantum rmmbcrs) . Finally, we observe 
that for fermion antisymmetry within the f3{uds) cluster, the /9(color-spin) functions must be symmetrized, as indicated in 
equation (7), since the flavor singlet carries the overall antisymmetry in \f3{u,d,s) >. 



B. Pentaquark Wavefunctions 

1. Color Octet (5(uds) States 

In this work, we restrict our attention to B = 1 molecules, which must have the flavor singlet 5(|^) as its baryon constituent. 



The color decuplet ^(|^) cannot form B = 1 molecules since there is no possible construction of color decuplet qqstates, which 
would be necessary to form physical, color singlets. However, we note that the C(|^) and its anti-particle ^^(|^) could form 
color singlet molecule states with baryon number 5 = (Cio — ?ro mesonic sextaquark) or _B = 2 (^lo — ^lo di-baryons). We 
write the explicit color-octet (spin-color symmetrized) flavor-singlet ^(5^) baryon wavefunctions labeled by color isospin {P,!^) 
and color hypercharge (Y"): |^c(spm^""'"")|/=,/|,y'= > 

Pt = \pt{\'')\l, ^, 1 > = 1^(8, l)|i^|rr6 > (12) 
= [uds — usd + sud — sdu + dsu — dus] (g) (13) 

^— [{2riribi — r-friftf — rir-fbf ) -|- Perms{l •<-> 3, 2 <-»• 3)] 
3v 2 

n1 = |nc(i^)|i, -i, 1 > = 1^(8, l)|^^|66r > (14) 
= [uds — usd + sud — sdu + dsu — dus] (g) (15) 

"^_ [(2b|b|r| — 6164. r| — 6|6|r|) -I- Perms{l ^ 3, 2 ^ 3)] 
3v 2 

S+ = |E+(1^)|1,1,0>= |^(8,l)|l^Hr> (16) 
= [uds — usd + sud — sdu + dsu — dus] ® (17) 

^^[(2rt5Tn - - nsrn) + Perms{l ^ 3, 2 ^ 3)] 

E» = 1^11(1^)11,0,0 >= m,l)\\^ \rhg{l) > (18) 
= [uds — usd + sud — sdu + dsu — dus] (g) (19) 

^^[(2rt6Tfl'i - n^ifft - n^'TSt) +Perms{l <-». 3,2 <-> 3)] 
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= |E-(i^)|l,-l,0>= |?(8,l)|i^|&s6> (20) 
= [uds — usd + sud — sdu + dsu — dus] (g) (21) 
^^[(2&|p|6| - b^g^b^ - 6j,5t&t) + Perms{l ^ 3, 2 ^ 3)] 

S° = |S°(i^)|i, ^, -1 > = 1^(8, l)|i^|ffffr- > (22) 

= [wds — usd + sud — sdu + dsu — dus\ (g) (23) 
v6 

^;^[(2ffT5Tn - S'T5i»'T " fliflT^t) + Perms{l ^ 3, 2 ^ 3)] 

H- = |Hr(i^)|i -1 -1>= |^(8,1)|1^|556> (24) 
= [uds — usd + sud — sdu + dsu — dus] ® (25) 
^;^[(2S'Tfl'T&i -STSi^T - 9i9]h) + Perms{l <-» 3,2 <-> 3)] 

A:! = |A«(i^)|0,0,0>= |e(8,l)|i^|r63(0)> (26) 
= [wds — usd + swd — sdw + dsu — dus] ® (27) 

^;^[(n^Tfl'i + n^'ifl'T + n^'Tfl'T) + Permsil <-> 2, 2 <-> 3, 1 <- 3)] 

(28) 

The color octet states are named according to their flavor octet counterparts. In this naming convention, the +, 0, — superscripts 
refer to the effective U{1) color charge of the state, not the electric charge (which is zero for any flavor singlet). To be explicit, 
we write the octet color electric charge (CEG) operator Q1 (for i = quark index): 

Qt = iU + \Y')i , (29) 

which is the color generalization of the U{1) electric charge quantization of SU{3) flavor octet states. The color isospin, 
hypercharge, and CEC quantum numbers are natural analogs of the corresponding SU{3) flavor quantum rmnibers. With 

the C/F transformation previously introduced, we have the following quark quantum number assignments for color isospin, 
hypercharge, and CEC as \color{I'',Ii,Y''\Q'') >: 

Miii||)> (30) 

IKi4i|4)> (31) 

|5(0,0,-||-i)>, (32) 

with the corresponding algebra: 

Q'\r> =l\r>, (33) 

Q'\b>=~\b>, (34) 

Q'\g> =-l\g>, etc. (35) 

Color confinement can be expressed as the condition that all physical hadrons must have total color electric 
chctrge equal to zero (and overall color antisymmetry). For all hadrons (including exotic multiquark systems) the total 
color electric charge is the sum over all quark contributions: Q'^ = Ql. 

2. Color Octet a ( us) States 

Next consider the a(gs)|8A, 2a|0~, 1^ > antisymmetric color-octet and flavor-doublet states {spin'"^'''*'")]!'^, I^jV^ >: 
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X+ = K+(0-)|l,l,l> (36) 

= -^{\rg > -\gr >) » -L(Ti - it) «> ^(9* - > (37) 

X° = |if°(0-)|i -i,l> (38) 

= -^ilbg > -\gb >) ® ^(Ti - it) ® - ^9) ' (39) 

7rc+ = k,+ (0-)|l,l,0> (40) 

" 71^'''^ ^ ® 71^^^ " -^^^ ® 71^''^ " ' ^^^^ 

7r° = |7r°(0-)|l,0,0> (42) 

= (^(l'-^^ > -1^^^ » - -^ilbb > -\bb >)^ (43) 
® ^(Ti-iT)®^(?s-sg), 

TTc" = kc"(0-)|l,-l,0> (44) 

" 71^''''' ^ "'""^ ® 71^^^ ~ ® 7!^^^~ ' ^^^^ 

^e° = |i^c°(0-)|i ^,-l> (46) 

= -^(|s6 > -\bg >) ® -^(Ti - it) ® ^(9« - ^9) . (47) 

K- = \K-{0-)\^,-^,-l> (48) 

= ^(lff^= > -1^=5 >) ® (ti - it) ® ^(9^ - «9) ' (49) 

J?" = l'?f (0^)10,0,0 > (50) 

= ^ ((|rf > -\fr >) + {\bb > -\bb >)) - 2i\gg > -\gg >)) (51) 

® ;^(U-iT)®^(9s-s?), 
where q = u is the positive electric charge I3 — and q = d is the negative electric charge Is = state. As with the color 



octet l3{B = 1) states, the color octet members of the a{B = 0) cluster are named in correspondence with their flavor octet 
counterparts. The spin one a-states are spin triplets having the same color and flavor structure \aT ispin'""''''')\r , li,^ >: 



K+ = \K+ii-)\^,^,i> 

= '^(\'^9 > -{gr >) <8) 

= |if:°(l-)|i -i 1> 
= -^{\bg > -\gb >) (» 



(TT),^(n + iT),(U) 



(TT),^(Ti + iT),ai) 



-j=(qs-sq), 



(52) 

(53) 

(54) 

(55) 



etc. for the other octet members. 
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3. Color Octet Pentaquark States 



With these Q;(qs) and /3(uds) basis states, we now form the totally antisymmetric (a — (3) molecules. We proceed by 
coupling color, flavor, spin, and space wavefunctions to overall quantum numbers of the physical (q — /3) states. As previously 
discussed, overall fermion antisymmetry in the /3-cluster follows directly from the antisymmetry of the flavor singlet F\{P). 
Construction of the symmetrized Fg{a(3\X) product function is straightforward yet tedious, so we omit the details. However, 
since antisymmetry of the pentaquark color function C\{aP\X) is crucial for color confinement, we show explicit details of the 
construction. To proceed, we note that each member of the anti-color triplet (f , b, g) is color-equivalent to the anti-symmetrized 
product of two color states: 

|f> =l=(\b>\g>-\g>\b>) (56) 
|6> =-^(|r>|s>-|5>|r>) (57) 
|5 > = ^ (|& > |r > -\r >\b>) . (58) 

When the mixed symmetry color function from the /3-cluster C£f g (/3) is coupled to the mixed symmetry color function from 

the a-cluster C^,jg{a), the SU{3) Clebsch-Gordon coefBcients for many pairings will be zero because those products can not 
contribute to a zero color charge state. For example, consider the couphng of a (rrb) color configuration in /3 with a (rg) state 
in a. If we symbolically express the color Clebsch-Gordon coefficient as CG(8a ® 8/3 -» Ix) = (a||Cf ,C'?||X^||Cf,C|,Cf ||/?) 
then 

Co {a\rg) <E) p\rrb)\X') = {a\\r,g\\X^\\r,r,b\\/3) = 0. (59) 

This coupling must be zero (and thus \rrb > (8)|rg > is prohibited) since the net color electric charge of such a pentaquark state 
is not zero: 

Q'' \P{rrb) (g) a{rg) > = {Qr + Qr + Qb + Qr + Qs) \l3{rrb) ® a{rg) > (60) 

= {1 + 1^1 + 1 + I) ^(^^'''''''^'^'^^''9) > (61) 
= (2) \l3{rrb) ® a{rg) > / 0. (62) 

Next consider a typical coupling between a — P color states that will contribute to a color neutral X-state: 

CG{a\rg)®0\ggb)\X^) = {a\\rrg\\X^\\g, gM\P) = ±^x, (63) 

where Nx is the X-state pentaquark color normalization, and the sign will depend on the adopted phase convention. Because the 
pentaquark color singlet wavefunction must have overall ^-symmetry, the normalization factor must have the same magnitude 
for all terms in the antisymmetric superposition. There are 9-primary (non-symmetrized) zero color combinations (only 8 
independent) with non-zero 8 (g) 8 — > 1 Clebsch-Gordon coefficients: {gf\rrb), {rg\ggb), (bg\ggr), {bf\rrg), {gb\bhr), {rb\bbg), 
{rf\rbg), (66|r6g), and {gg\rbg). Explicit color antisymmetry can be constructed in a sequence of color transformations starting 
with an arbitrary primary pentaquark color combination. For example, starting with the (ai=i|/9i=i) = {gf\rrb) term, we first 
write the product of (non-normalized) MS-|a > and MS-|/3 > color octet functions: 

Clts(Mgr\X) Clts{Pi\rrb\X) = (,gf - fg) {rrb - rbr + brr) , (64) 

and hkewise for the 8-independent products. Next transform this first term by applying the color transformation r <-> & (and 
subtract to anti-symmetrize): 

Clis{gh\X) Clts(bbr\X) = (gf - fg) {bbr - brb + rbb) . (65) 

Repeated application of the various color transformations: r<r-*b, r^g,b<--*g generates the following (normalized) antisym- 
metric color cluster product expansion: 

8 

C\{al3\X) = M^x ^(-1)'+' Clts{ai\X) Cltsmx) , (66) 

i=l 

where the la >-color functions are: 
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Clis{ai\9r\X) = {9f-fg) (67) 

Clts{cx2\gh\X) = {gh-hg) (68) 

Clts{az\rh\X) = {rh-lyr) (69) 

Clts{ai\hf\X) = {hr-fh) (70) 

Clts{ccr,\r-g\X) = {r-g~-gr) (71) 

Clts{m^bg\X) = {hg-gh) (72) 

Cl(a7|rf, 66,flg|X) = {rf - fr + bb -bb + gg - gg) (73) 

C|(a8|rf, &6, ffffIX) = (rf + fr + 66 + 66 + 55 + 35) (74) 



and the \/3 >-color functions are: 



CMsiPi\rrb\X) = {rrb - rbr + brr) (75) 

Clis{l32\bbr\X) = {bbr - brb + rbb) (76) 

Cms (A 1 665 1 X) = {bbg - bgb + gbb) (77) 

Chs{l34\rrg\X) = {rrg - rgr + grr) (78) 

Cltsm99h\X) = (ggb - gbg + bgg) (79) 

Cl,sm99r\X) = {ggr - grg + rgg) (80) 

C%{(37\rbg\X) = {rbg + rgb + grb + gbr + bgr + brg) (81) 

C\{l3&\rbg\X) = {rbg - rgb + grb - gbr + bgr - brg) (82) 

Note within both the a and l3 octets there are two terms that have zero net color charge. These color neutral octet member 
could be written with MS-symmetry, but to make the Clebsch-Gordon factors simple (and to demonstrate explicit A-color 
symmetry) we show the two independent S and A-symmetry combinations. After the sum in Eq.(66) over the 8-octet states is 
performed, the color function will be totally antisymmetric in color interchange. This means that the interchange of any two 
colors will generate the opposite overall sign in the color cluster product expansion. For example, under the color exchange 
r •<-» 6: 

Cl,s{ai\gf\X) C^Ms{0i\rrb\X) ^ (-)C^s(a2|ff6|X) Ct,s{02\bbr\X) (83) 

Ct,s{a3\rb\X) Chs{3s\bbg\X) ^ {~)Ciis{a4\br\X) Ciis{04\rrg\X) (84) 
Clis{ao\rg\X) Ciis{l3o\99b\X) ^ {-)Clis{a,\bg\X) CtiAW6\ggr\X) (85) 
cUa7\rf,bb,gg\X) cUP7\rbg\X) ^ {-)C%{a7\bb,rf,gg\X) cUP7\brg\X) (86) 
CUa8\rf,bb,gg\X) C%{^8\rbg\X) « (-)CUa8\bb,rf,gg\X) C%{^8\brg\X) (87) 



The QMM pentaquark color antisymmetry is now explicit. There are a total of (6 x 2 x 3) + (2 x 6 x 6) = 108 terms in the 

_L _ 1 

/108 ~ eVs' 



antisymmetric sum, hence the color function normalization is: Ax = ^ ~ ^ 



C. Color Magnetic Moments and Confinement 

We now turn our attention to a dynamical constraint that has two novel consequences: (1) to provide a selection rule that 
limits the QMM excitation spectrum (allowing only 4-physical states); (2) to provide dynamical relations between constants of 
the QMM Hamiltonian such that our calculated pentaquark baryon mass spectrum has only one free parameter (the ^ mass). 
Using the color electric charge ( GEC) operator described in the previous subsection, we now define the color magnetic moment 
(CMM) operator: 

A" = ^ A< o-3i = /"o ^ Qi *3i , (88) 

i i 

which is the color generalization of the U(\) magnetic moment operator of the conventional (non-exotic) quark model. The 
sum extends over all quarks and the color magnetic moment operator /t? satisfies the SU (3) algebra: 

jiy > =nlQ''\r>= I |r > = Mr |r > , (89) 

fi^\b > =nlQ^\b>= -'^^j.l\b>= ^l\b> , (90) 



A'^ls > = Q'^ Is > = - o Iff > = Ms Is > • (91) 
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The constant /iq is the fundamental color magnetic moment, which may be regarded as a free parameter of the QMM, however 
we take a different approach since its numerical value is not used in our spectrum calculations. The important point about the 
constant fig is that its value (proportional to the strong coupling constant Qs = \/^nal) must be the same for all quark colors 
to satisfy exact SU (3) color symmetry. 

We can now state our hypothesis of color magnetic confinement (CMC): 

The net color magnetic moment is zero for all physical hadrons. 

If the CMC hypothesis is true, then the converse must also be true: 

Any hadron with non-zero color magnetic moment is unphysical. 

To see why the CMC conjecture is reasonable, we apply our hypothesis to a conventional quark model flavor octet state (the 
proton). First we take the direct product of standard A-color and MS-spm. functions for the proton: 

|P :T>2-3 = -^{rbg - rgb + bgr - brg + grb- gbr) (92) 

X ^ T (Ti - it) , (93) 

where the (2-3) subscript denotes A-symmetry for interchange between the second and third quarks in the spin function (and 

no other pairwise interchange symmetry). Obviously, the proton has CEC equal to zero. After applying the CMM operator, 
there are a total of 36 terms (12 terms with 3 quarks each) with exact cancellation between terms: 

<T: p\ ii" \p :T> = <T: p\ ^ Ql ^si \p :T> = . (94) 

i 

Note the cancellation is only exact if /Uq is the same for each color. The physical proton is built from a superposition of 

MS'-symmetry flavor functions, but the exact cancellation in the CMM holds separately for each one, so the proton satisfies the 
CMC hypothesis. We can show that every member of the (non-exotic) baryon (^^) octet, (|^) decuplet and meson (0^, 1^) 
nonets also satisfy the CMC hypothesis (but omit the details here). It is possible that some (or many) of the excited states of the 
conventional (non-exotic) quark model may not satisfy the generalized CMC sum rule (including orbital CMM contributions) 
and therefore drop out of the excitation spectrum. The CMC principle may help to resolve the longstanding missing nucleon 
(A'^*) resonance problem of baryon spectroscopy, or CMC may simply not be generally true, however we leave this question 
open for future consideration. We now proceed to apply the CMC conjecture to constrain the Hamiltonian of our pentaquark 
baryon molecule states. 



D. Color Magnetic Confinement (CMC) Constraint 

The complication of applying CMC in the QMM pentaquark basis is that there must be at least one unit of relative orbital 
angular momentum between color octet \a{qs) > and \[3{uds) > clusters. Hence, there should be a CMM contribution from 
the colorful orbital motion of the molecular state. If CMC is generally true, then we must calculate the orbital contribution to 
the total molecular color magnetic moment. When this is accomplished, there will be consistency relations between the color 

spin-spin and spin-orbit coupling strengths in our QMM Hamiltonian so that overall CMC is satisfied. 
Our molecular pentaquark spectrum is generated by the following Hamiltonian: 




where S — Sa + 8/3 is the total spin, and L is the orbital angular momentum operator. This Hamiltonian has the same form as 
a non-relativistic reduction of the Dirac equation with a minimally coupled gauge field 4-vector potential (generating electric 
and magnetic interactions). Cluster confinement is enforced by the harmonic oscillator potential, which generates a convenient 
S.H.O. basis of radial wavcfunctions for our spectrum calculations. As usual, is the strong coupling constant of QCD [14] 
(fixed by conventional quark model spectroscopy as = 0.3). We checked that using the running coupling as{r) gives an effect 
less than 1% in the radial expectation values calculated in this work. The factor C% multiplying the Coulomb term is the 
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color factor of one gluon exchange between color octet quark clusters [14] (derived later). The Aa/s and Bap coefficients are 
combinations of color gyrornagnetic ratio 7-factors (analogous to the definition in electrodynamics) corresponding to the color 
magnetic moments of the quark clusters and orbital motion, which are constrained and related by the CMC sum rule. 

We proceed by defining the total color magnetic moment of the molecular pentaquaxk system jltot in terms of the intrinsic 
fla and fig cluster contributions, and the orbital contribution jTj^: 

fitot = fia + + ■ (99) 

The generalized CMC constraint is: 



CMC: Pl.t = 



The intrinsic a and (3 magnetic moments are expressed: 



where gs ~ V Airas is the QCD strong coupling constant, and the (70,7/3) factors are dimensionless color gyrornagnetic ratio 
factors of the (q,/3) clusters. It is possible (yet tedious) to relate the color gyrornagnetic ratio constants to the fundamental 
color magnetic moment /Ltg of the CMM quark operator using the explicit cluster basis wavefunctions of section II B (analogous 
to how traditional quark model SU{3)f baryon magnetic moments are related to quark magnetic moments). However, we 
apply the generalized CMC rule to constrain the effective molecular CMM constants, and hence do not pursue these QMM 
parameter relations in this work. Note the baryon /3 cluster does not have a factor of 2 in the mass denominator because we 
assume the Landau g = 2 factor of Dirac electrodynamics which enhances fermion magnetic fields relative to non-fermion fields 
(such as those generated by \a >). Next, we define the intrinsic (a/3)-molecule CMM and its corresponding effective color 
gyromagnetic ratio 7^'^ by coupling to total spin: S = So, + Sg: 

■.J-S a/3 



Ma/3 = Ma + = i-y-' 7r ( S , (102) 

where Mag = ^"^^^ is the (a-/3) reduced mass, and the J-dependent phase is required for a technical reason we come to 

shortly. The orbital color magnetic moment /tt^ is expressed in terms of the total orbital angular momentum L and effective 
orbital gyromagnetic ratio 7^'': 

The origin of the orbital gyromagnetic ratio factor can be traced to the fact that the total orbital color magnetic moment is 
generated asymmetrically between the colorful clusters of different mass. Each cluster carries a different amount of orbital 
angular momentum relative to the center of mass of the molecular system (Fig. 1). 

For orbiting (a-/9) clusters in a particular spin configuration, each carries individual angular momentum: 

f ^ ^ classical 2 \ ^ 

La = Ta X Pa [Ma Ta LOl) Z (104) 

=> ^ ^ classical 2 \ /^^^n 

Lg=rgX pg — (M/3 rg UJl) z (105) 

where z is the quantization axis of spin and orbital angular momentum, and the angular velocity (expectation value) coi must 
be the same for both clusters so the center of mass position is fixed in space (i.e., to satisfy conservation of linear momentum 
with no external forces). Also, the quantum two-body system satisfies Newton's 3^'' Law (with expectation values used for 
radial factors): 

Ma Ta Uj'i = Mg rg UJl . (106) 

Total orbital angular momentum is the vector sum of both cluster contributions, and naturally generates the standard reduced 
mass Mag factor of a CM separated two body system: 

L = La + Lg (107) 
= {Ma rl UJl + Mg rg uji) z (108) 
= Mag (r„ + rgf uji z . (109) 
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Because of Newton's 3'^'' Law, Eq. (106), each cluster angular momentum is proportional to the total angular momentum, 
depending on mass ratio factors: 

^=©4 (110) 
^« = [1 + (M./M,)] ^ 



Mo 

M,}) [1 + {Mo^/Mg)] 



(112) 



Returning to the orbital color gyromagnetic ratio, we couple the separate a and /3 orbital CMM contributions to the total 
orbital angular momentum: 



V 2M, 



Using Eqs. (111,112) we obtain: 

"^^ (M„ + M^)2 



Mg Ma 



(114) 



It is interesting to note the equal mass and infinite asymmetric mass limits: 



lim. 7f-i (115) 



Ma, = M0 2 



ct0 



lim 7^ ^1 (116) 



lim 7r ^ 1 • (117) 

We now explore the effective color gyromagnetic ratio of the coupled molecular spin CMM, 7^^^. Using Eqs. (100,101,102), 
we have: 

The J-dcpcndcnt phase factor is required for consistency because in the coupled quantum mechanical molecular system, 
(Sa, S/j, S, 5*2, L, Lz, J, Jz) arc "good" quantum numbers (operators commute with Hamiltonian) , but {Saz,Si3z) are not. To 
obtain the required relations, first examine the classical +2 projection of Eq. (118) with {S = ^\Sa = 0, Sg = |): 

Our first consistency relation for 'yg'^{S\Sa, Sg) follows directly: 

i-y-^yfiW.) = (^) 7.., (120) 



' 2 ' ' 2 ' \ M^^ ^ 

where a = Kc{qs\0^) (a spin 0^ color octet cluster with kaon spin-flavor composition), and f3 = ^c{uds\^^), the previously 
discussed color octet, flavor singlet baryon cluster state. Because the Kc{0~) cluster has spin zero, its color gyromagnetic ratio 
must equal zero: fK^ = 0. 

Next consider a spin one (color octet) meson cluster {Sa = 1), the a = K*{qs\l~) with vector kaon spin-flavor composition 
coupled with the /3 = ^ciuds\^~^): 

f (ill, i) = M,„.(^-^), ,m) 

The last possible molecular spin combination is the vector a = K*{qs\l~) coupled to the /3 = ^ciuds\^^) in a S = | 
configuration: 

^J-i „c«/3.3,, 1, _ ,1, / 7ifJ , 7€c 



i-y-- 7^(^11, o) = (3) ^-cCc (W^ + ^J- (122) 



11 



Finally, we apply the CMC constraint to develop eonsistency relations between the effective gyromagnetic ratios of the cluster 

states (7ifc,7if* ,75) and the orbital factors. Eq. (99) can bo written as: 

= + = 7f (^) S + 7f (^) £ (123) 

f gs \ f CMC 

So the CMC condition requires the color gyromagnetic ratio of any physical molecular state with quantum numbers 
{J,Jz\L,Lz,S,Sz) to be identically zero: 'y^f^i^s) ~ ^- Examining the z projections for the different {J\LS) quantum cou- 
plings, we find the following consistency relations for al3{J\LS): 

Kci{\\l,\): 7f(i|0,i)=7^« (125) 

^4(f|l,i): -7f(i|0,i) = -7f=« (126) 

KU{\\1,\): 7f (^|l,i)=7f*^ (127) 

-'if{\\h\) = -lT^' (128) 

^c^(i|l,^): -7f(^|l,^) = (i)7^*' (129) 

^:e(f|l,f): 7f(f|l,f) = (130) 

^c*e(f|l,f): -7f(f|l,f) = -(i)75' (131) 

Note the only consistent solution for Eqs. (129,130,131) is: 79'^(||1, |) = 0. But that would also require 7^'°^ 0. Since 

7^"''^ is uniquely determined by Eq. (114) (non-zero), the K*£^{J\LS) — (||1, |), (||1, |) states cannot consistently satisfy the 
CMC condition and are therefore unbound (unphysical states) in the QMM. Careful inspection will show the sign difference 
between the J = ^ and J = | configurations of the L = 1 and S = ^ Kc£, and solutions correspond to the J-dependent 
phase factors in Eqs. (120,121). Now using Eqs. (120,121,122), together with the two independent (CMC consistent) solutions 
Eqs. (125,127) (or Eqs. (126,128)) we find the corresponding set of relations between the QMM gyromagnetic ratio factors: 

7?(||0.|)-7f-'- = (^) 7t. (132) 



which we invert to express the unknown (phenomenological) cluster color gyromagnetic ratio factors 7if j and 75^. in terms of 
the derived orbital factors 7™'^ , and colorful cluster masses Ma = {Mk^ , Mk* ) , Mjs = M^^ : 

It is important to stress here that these parameter relations were derived analyzing the cleissical z projection of the coupled 
spin Sz = Saz + Spz, but we have already noted that Saz and Spz are not good quantum numbers of the molecular system. 
We will refer to this parameter set as the "classical approximation" for the QMM-CMC constrained gyromagnetic ratios. 

However, we can derive a set of parameters consistent with the rules of quantum mechanics simply by taking expectation 
values of the spin operators cissociated with Eq. (118). Using the standard spin (and later orbital) expectation values: 

<S-S> = 5'(S'+1) (136) 
< S; • 51 > = SaiSa + 1) (137) 
<Si3-Si3>= S0{Si3 + l) (138) 

<S'a-S0>= i[5(5 + 1) - 5„(S„ + 1) - SgiSg + 1)] (139) 

<L-L>= L(L + 1) (140) 
<J-J> = J(J+1) (141) 

<£■§> = i[J(J + 1) - L{L + 1) - SiS + 1)] , (142) 
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by taking the S • S projection of Eq. (118), we obtain the general quantum consistent parameter relations: 



(143) 



S{S + 1) 



s.is. + 1)(^) + s,is, + 1)(^)+ < -s.Xj^^+if;, 



Applying Eq. (143) to the cluster a = Kc{0 ) {Sa = 0, 5/3 = i, S = 1) (and canceling the phase factor when coupling to 

2' 2 



different total angular momentum = | , | states) , we obtain: 



M, 



(c 



7€c 



(144) 



which is the same result as the classical approximation Eq. (132). A difference shows up when we apply Eq.(143) to the cluster 

a = K*{1~) {Sc, = 1, Sf3 = 2, S = I): 



a/3/ 1,-1 Ix^-fJSc _ ( 2 7k; . , 



1 7?c 
3M5, 



^2' '2' 

Inverting as before, we obtain the quantum and CMC consistent gyromagnetic ratio parameter relations: 



7So = 



7k, 
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2 

7Kc = . 



7i 



2 Mk, 



7i.' 



(145) 



(146) 

(147) 
(148) 



E. Color Electro-Dynamics (CED) QMM Hamiltonian 

We now relate the color gyromagnetic ratios to the Aag and Bag coefficients in our QMM Hamiltonian Eq. (95). In the 
spirit of our CEC and CMC conditions, we treat the form of the color octet cluster spatial and spin dependent interactions 
assuming that the primary non-Abelian features of QCD are contained in the confinement term, leaving the spin-spin and 
spin-orbit terms in the same form as the usual electrodynamics electric and magnetic interactions. First consider the spin-orbit 
interaction of electrodynamics for a charged boson a and fermion /3. Each cluster separately contributes to the total spin orbit 
interaction: 

H^^ = HZ + Hi (149) 
.1^ 1 /J^\ fi^\ 



As already noted, the magnetic moment of a fermion is twice as large as a corresponding boson. The factor of 1 in each term is 
due to the Thomas precession effect. The spin-orbit interaction is derived from the energy of a magnetic dipole in a magnetic 
field B: 

H,„ = -jl-B (152) 
P=ij^)S (153) 



B : 



^MaB \r dr J ' 



Here g is the fundamental charge (electric or color), Alag is the reduced mass of the two body system, and for a Coulomb 
potential the derivative term generates a radial dependence. We write the derivative term explicitly so that the QMM 
color generalization for the radial dependence can be modified by the confining term in the potential Eq. (97). Using S.H.O., 
confinement generates a constant term (the classical "spring constant") in Eq. (98), but the extra term would involve a 
term for exact linear confinement. We define the effective spin-orbit interaction in terms of a coefficient Bap: 



B. 



a/3 



S = 



Ma Ma \ 2M, 



(155) 
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then the total spin-orbit interaction takes the form expressed in the QMM Hamiltonian Eq. (95). Matching coefficients, we 
define two constants (a, 6) which depend on the gyromagnctic ratios (70,7/3) and mass ratios: 



Baf3 S = a Sa 



bSp 



V 2M^ J \ 2 J 



b = 



2M, 
2M« 



7/3 



(156) 
(157) 

(158) 



Projecting both sides of Eq. (156) with the total spin operator S, and evaluating the spin expectation values of Eqs. (136-139), 
we find the spin-orbit Hamiltonian coefficient: 



{a + b) + {a- b) 



S{S + 1) 



(159) 



Note the Bap coefficient depends on the spin quantum numbers {S, Sa, Sp) of the molecular system. 

Next consider the more complicated spin-spin (SS) interaction. The form of the SS interaction is different for the L = 
and L 7^ cases. For completeness, we consider both terms. Our QMM is restricted to L = odd, but we investigate the L = 
Hamiltonian to see how it aifects our dynamically restricted spectrum. The electrodynamics inspired form of the interaction is: 



-i--- (?) (^ 



7/3 
Mg 



ldV\ 
r dr J 



(160) 
(161) 



In the case of L = 0, the expectation value of < r ' > in the derivative term diverges. However, the spin-spin term gets a 
(5-function correction which regulates the divergence, hence the < >l=o factor is replaced by [14]: 



Roo{r = 0) 



(162) 



where Rooir) is the normalized radial wavefunction (in the n = 0,L = state) described in the next section. In the case of 
L ^ 0, we must take the expectation value of the spin operator in Eq. (162). To simplify the evaluation (and recover the form 
of the QMM Hamiltonian) we note that to preserve the CMC constraint with non- variable gyromagnetic ratio parameters, the 
spin and orbital angular momentum must couple in maximally aligned or aiiti-aligned configurations (otherwise, the dipolc 
color magnetic fields of the quark cluster will not exactly cancel the orbital dipole field). Because the CMC condition locks the 
physical color singlet molecule into maximally aligned (anti-aligned) configurations, all quantum substates with zero angular 
momentum projections (i.e.. Ml = O.Ms,^ = 0,Ms^^ ~ 0,Ms = 0) are prohibited in our QMM. This physical restriction 
(required for consistency with CMC) greatly simplifies the spin-spin operator, since the expectation value of any Sx or Sy spin 
projection will be identically zero: 



[3(5„ • r){Sp -r) - S^-Sp] 2 5„ • S^j . 

So we arrive at the QMM spin-spin Hamiltonian term expressed in Eq. (95). The Aq 
gyromagnetic ratios for the a and /3 quark clusters: 



(163) 

is therefore a simple product of the color 



Aal3 = Jail} ■ 

We now discuss the QMM spatial wavefunctions and resulting mass formula for the molecular pentaquark exotics. 



(164) 



F. QMM Spatial Wavefunctions and Expectation Values 

To determine the molecular exotic baryon spectrum, we calculate the first order energy corrections to the exact energy 
eigenvalues associated with the confining harmonic oscillator Hamiltonian Hq: 

where in this subsection /i is the reduced mass (/x — Map = MaMp/{Ma + Mp)), and w is the quantized oscillator frequency. 
In this work, we use natural units where Ti = c = 1. The angular eigcrifunctions of Hn arc the spherical harmonics, Yira{d,(f))- 
The energy eigenvalues of the unperturbed oscillator are obtained from the radial Schrodinger equation which includes the 
centrifugal barrier term: 
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X + 

where 



X = 0, (166) 



k = y^2^E , and X = ^lu , (167) 

x(r) = rR(r) is the reduced radial wave function. Following the usual procedure, the small/large-r asymptotic behavior of 
Eq. (166) determines the form of the general solution: 

X(r) = r'+^e-^^'w(r) , (168) 
and the function u(r) satisfies the differential equation: 

i u"{x) + + + EHu{x) = , (169) 

expressed in terms of the dimensionless radial variable x = Xr^. The regular solution of Eq. (169) is a confluent (or degenerate) 
hyper-geometric function ^[a,b;x] [15,16]: 

$fa 6- = V + x' (170) 

na, *J 2^ r(a) r(6 + j) r(i + j) ^ 

j=0 

_i , «^ , aia + i) x^ , a{a + l){a + 2) x'' , 

~ ^ 6^ 6(6 + 1) 2! ^ b(6+l)(6 + 2) 3! ^ - ^ ^ 

«(^)='J>[^(' + ^-^M, ^+1 (172) 

Since physical solutions must have normalizable wave functions, the infinite series in Eq. (171) must be truncated after a finite 
number of terms a = —n (integer n), and hence the oscillator energy levels are quantized: 

l(l+^-E/co) = -n (n = 0,l,2,3,...) ■ (173) 

Solving Eq. (173) we find the harmonic oscillator energy eigenvalues: 

Er,i = {2n + l + ^)u; , (174) 

associated with the eigenfunctions: 

•fnim = N„, r' e-T-"' $[-n, / + ^; Ar'] Yi^{9, 4>) ■ (175) 

With these basis functions, we compute the expectation values of the full Hamiltonian to determine the molecular exotic mass 
spectrum. In this work, we consider only the states with n = 0. Orbital excitations are prohibited by the CMC constraint, 
however radial excitation may bo permissible. Wc do not explore the radial excitations and their consistency with CMC in this 
work. Since $[0, 6;x] = 1, the normalized wave functions are: 

(1+2 I — \ ^ 
(2/^+1)!! ^'^ V^J (1^6) 

Using these normalized wave functions, a general expression for the expectation value of the radial coordinate raised to an 
arbitrary power [p) can be evaluated: 

/poo 
dn J r^+'' dr \^oimf (177) 

(3+p) 



(2/ + 1)!! Ytt 



Specific expectation values for / = 0, 1 and p = 2, 1, —1, —2, —3 are evaluated in Table I. 

The mass formula for the exotic molecular systems associated with the Hamiltonian of Eq. (95) is therefore: 
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Mf =M„ + Mg + {2n + l + ^)ij + C%a, < r"' > (179) 



A. 



al3 



MaMp 



<Sc,-S(3> [-C% <r-^> + ] 



with Sa = 0, 1 and sp = |, |. The oscillator frequency uj in Eq. (179) is given by the mass sphtting between the proton and its 
first radial excitation using Eq. (175): 

w = ^[M;v(i440) - Mp(938)] = 0.251 ; (180) 

The usual color factor between color triplet quarks forming a color singlet Cf(3 (8> 3 ^ 1) = —4/3 is not appropriate for 
onc-gluon exchange between color octet clusters. Therefore, we will derive the correct factor Cf. — Cf(8 x 8 ^ 1) for octet 
bonded clusters. Evaluating the Feynman diagram for gluon exchange between color octet (point-like) particles (Fig. 2), we 
have: 

C% = i^f (i)^ n'" f^'^' 5,,, 5^f, S^,p, (181) 
= -(l/8)(6)[ {r'f + 6 {f'^'f + 2 {f^y ] ^ -3 , (182) 

which involves a trace over the product of SU{3) structure constants f^i^^^ and projecting octet color functions (ai^) onto 
singlet initial |1 > and final |1' > state color functions: 

|1 > = (aia^)^ = -^^c.ff (183) 

|1'>= {al,4,)^ = ^ W ■ (184) 

Note the color octet-octet interaction generates a strongly attractive one gluon exchange potential, Vig{r) — C%as/r. 

The last step in constructing the QMM spatial wavcfunctions involves coupling cluster spins (Sa,S[j), total spin S, and 
orbital angular momentum L to overall angular momentum J. The angular momentum J-states are calculated using the 
superposition of |L > |iS > basis states with CMC restricted Clebsch-Gordon coupling factors (all magnetic substates equal to 
zero are prohibited in the CMC restricted sum): 

CMC 

\'i^,fi{r,e,4>)\\S,L,ML,J,Mj >= ^ (185) 

C{J,Mj\L,Ml,S,Ms) C{S,Ms\Sc,,Ms^,S0,Msp) RNL{r) Y^^{e,4>) . 

Since the CMC principle restricts the angular momentum substates, the standard Clebsch-Gordon factors are renormalized 
in the QMM sum over states (details omitted). We now discuss our QMM parameter assignments and the resulting physical 
spectrum of pcntaquark molecules. 



G. QMM Parameters and Spectrum Results 

The QMM cluster masses are not a priori known (since the color octet clusters can not be measured in unbound isolation). 
However, their masses should be comparable to the physical K, K* and A masses due to the dominating contributions from 

constituent quark masses. Since the quarks inside each cluster combine to form a color octet, we anticipate a mass shift (increase) 
originating from the repulsive color factor of one gluon exchange between quarks in the octet color configuration. To estimate 
this mass shift, we calculate the color factor of one gluon exchange between qq coupled to a color octet, Cf(3 (8> 3 ^ 8) = +|, 
whereas the color factor for the color singlet is Cf(3 (g) 3 ^ 1) = — |. So the estimated mass shift due to one gluon exchange 
between quarks inside a color octet cluster is: 5M = |as < 7 >. 

For a typical hadronic size (< r >« 0.5 fm), we deduce the mass shift 5M ~ 0.1 GeV. Since this is a small number compared 
to the bare cluster masses (i.e., physical K,K* and A masses), and because the ^-mass is an unknown adjustable constant 
which is fixed by the lightest pcntacjuark candidate, we simply use the physical K, K* masses for the a{qs) clusters aud absorb 
the mass shift ambiguity into the unknown phenomenological ^-mass, which is the only unconstrained parameter of QMM 
(Table II). The calculated QMM masses of the exotic (pentaquark) states are summarized in Table III. 
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III. EXPERIMENTAL STATUS 



In this section, we review the current status on experimental data focusing on associated electromagnetic production of 
strangeness systems of protons. In particular, (re-) analysis of experiments were done in the spirit of finding potential narrow 
structures to support our QMM calculation. 

A. Photo-production channel 

Guided by the recent claim of three narrow resonances from the SPHINX collaboration [18-20] , we re-examined the published 
data on K-meson photo-production using the p{-y, K)[A, E°] reaction from the SAPHIR detector at ELSA [21]. 

The invariant mass distributions for both A and E hyperons are compared to the background generated by a well constrained 
hadrodynamical Lagrangian model prediction of [2,3] in Fig. 3. The various curves account for when the narrow baryon X- 
resonances are included (solid line), or not included (dot-dashed line). We assume assume relativistic Breit-Wigner shapes for 
the resonances, which are superimposed on the smooth background. The masses, total widths and signal strengths are adjusted 
to optimize agreement with the analyzed data (Table IV). 

Overall, the qualitative agreement between the data and the Lagrangian model prediction is significantly improved for both 
hyperon production channels when the X-baryon resonances are included. A detailed description of the hadronic Lagrangian 
model employed is under preparation for submission [2,3]. 

B. Electro-production channel 

The kaon photo-production results described in the previous section provided motivation to investigate the equivalent (virtual 
photon mediated) kaon electro-production reactions. The elementary process for p{e,e' K)[A,'SP] is well described at low 
[22]. Analysis of recent Jefferson Lab data has shown that isobar based models provide a good description of electro-production 
observables [23,24]. We employ the same isobar model [2,3] to check for consistency with the photo-production results. We 
analyze dedicated kaon electro-production data that were taken in December 1999 during the second stage of experiment 
E91016 [25,26] at Jefferson Lab in the experimental Hall C using the 100% duty cycle of the CEBAF machine at a fixed beam 
energy of 3.245 GeV. 

The scattered electrons were detected in the High Momentum spectrometer (HMS, AP/P = ±10%, 5P/P = 10"^, AQ = 
6.7 msr), set at 19.64° and with a central momentum of 1.432 GeV/c. The K-mesons were detected in the Short Orbit 
spectrometer (SOS, AP/P = ±20%, SP/P = 10"'\ AO, = 7.7 msr) along the virtual photon direction at 14.23° with a 
central momentum of 1.290 GeV/c. The corresponding invariant mass and settings were 1.925 GeV and 0.55 (GeV/c)^, 
respectively. The detector packages in both spectrometers are similar and include: two planes of drift chambers (used for 
tracking information), four planes of plastic scintillators (used for trigger formation, and time-of-flight (TOF) information), a 
gas Cerenkov (for e/ir separation), and a lead glass shower counter (for electron identification). The SOS had two additional 
Ccrenkov counters to improve identification of the kaons: an aerogel Cerenkov (for n'^ /K'^ separation), and a Incite Cc^reukov 
(for p/K^ separation). The detector performances achieved were similar to the ones obtained in earlier kaon experiments 
performed in this experimental hall [23,24]. 

In the hadron arm, the K'^ were identified using the aerogel, the Incite and TOF. In the electron arm, we used the gas Ccrenkov 
and Pb-glass. A sample of 8262 coincident kaons were then reconstructed which includes subtraction of the background events 
from the target cell walls (0.6%) and random coincident events (rcal-to-random ratio about 10/1). 

The two hyperons produced during the reaction can be seen at their expected location (Fig. 4): Ma(iii6) = (1116.40 ± 
0.66) MeV [Fa = (4.82 ± 0.65) MeV] and Ms0(ii93) = (1192.7 ± 0.12) MeV [Fe = (4.86 ± 1.60) MeV]. 

The invariant mass distributions are shown in Fig. 5 in addition to our model differential cross section of [2,3] for the A (left) 
and (right) hyperons. The fits of the three enhancements using relativistic Breit-Wigner shapes are listed in Table V. 

As a result of the different bin sizes used in the analysis of the data in [21] (30 MeV) and [25,26] (~12.55 MeV), the new 
Jefferson Lab electro-production invariant mass distributions axe able to resolve two nearly degenerate resonances, X\ and X2, 
in the vicinity of the SAPHIR X2 [21]. We believe that due to the large bining of the SAPHIR data, the first peak in E" 
production at 1.974 GcV averages over the JLab Xi and X->. which contribute in the same SAPHIR energy bin. The resulting 
averaged value, = (1974.1 ± 6.7) MeV [F ^^^^^^ = (97.7 ± 55.1) MeV], is less than 1-cr of the extracted SAPHIR 

X2 mass value, and the decay width is approximately equal to the sum of the two JLab Xi and X2 resonance widths. As 
in photo-production, we find optimum agreement between our model prediction and the electro-production data when the 
X-resonance degrees of freedom are included [2,3]. 
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C. Summciry 



Although the resonances extracted have fairly poor statistics, there is a clear definite pattern between all of the experimental 

data. This gives us a strong confidcnco for the existence of tlu' narrow states extracted. We summarize the results or our 
analysis for all experimental data sets in comparison with the QMM spectrum predictions in Table VI. 

Although a dedicated new high precision electroproduction experiment is planned to be performed for confirming and 
extracting the properties of these resonances [28], additional data in photo- and electro-production in that W range will be 
needed in the near future to shed light on the possible existence of the combined SPHINX, SAPHIR and Hall C data resonances 
discussed here. We would like to stress that the W distribution of the differential cross section would have to be kinematically 
constraint to center-of-mass angles between the virtual photon and outgoing kaon to be < 20 — 30° [2,3]) due to the 

(expected) strong forward peaking effect of these states on observable quantities. 

IV. CONCLUSION 

We have developed a quark molecular model (QMM) of S = pentaquark baryons assuming color-octet bonded (qs)-(uds) 
molecular structure, and predict the low lying mass spectrum (4-states) using a color electrodynamics (CED) Hamiltouian. 
We have introduced the idea of color magnetic confinement (CMC) analogous to the well known color electric confinement 
(CEC) of traditional quark models. The CMC condition limits the allowed physical states in the QMM spectrum and generates 
consistency relations between the color gyromagnetic ratio factors of the QMM-CED Hamiltonian. This novel prediction stands 
in contrast to the usual outcome of quark model spectroscopy (which generate an overabundance of excited states) and provides 
a possible explanation of the dilemma of '^missing resonances" . 

To show the plausibility of our pentaquark spectrum predictions, we have reported experimental evidence of four narrow 
structures in the vicinity of = 2 GeV in associated strangeness production reactions using the recent ELSA data in kaon 
photo-production [21] and Jefferson Lab data in kaon electro-production [25,26]. One resonance structure is not consistent with 
the QMM mass prediction and requires a typical (large F ~ 200 MeV) hadronic decay width, and hence is assumed to be a 
conventional (possibly missing) A'^* resonance, the A^*(1898). 

The suggested candidate pentaquark states of this analysis have average masses of 1.812, 1.962 ^, 2002 and 2.064 GeV, all 
of which have been identified with anomalously narrow decay widths (few tens of MeV). Two of those resonances (1.943 GeV 
and 2001.7 GeV) seen in the recent JLab data are believed to appear as an unresolved single resonance in the SAPHIR data 
at 1.974 GeV (due to relatively coarse energy bining). The fitted Breit-Wigner resonances align closely with the 4-predicted 
pentaquark states of the QMM, with masses: Xo(1801)[l/2+], Xi(1967)[3/2+], X2(1999)[l/2+], and X3 (2064) [3/2+]. 

Additional high resolution experiments are required to establish the unambiguous existence and properties (i.e., masses, 
decay widths, branching fractions, hadronic coupling strengths, form factors, etc.) of the candidate 5 = pentaquark baryons 
suggested in this work. Our results provide a strong indication of possible crypto-exotic pentaquark baryon resonances in 
associated strangeness production channels, and support the need for investigating reaction channels other than non-strange 
(tt, p, w, etc.) production to look for missing quark model resonances [29] and QCD allowed exotic states. 
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TABLE I. Expectation values of the radial coordinate in the oscillator basis. 



Cluster {a, (3) 


Cluster Mass 


Gyromagnetic Ratio 7 


Gyromagnetic Ratio 7 




[GeV] 


(Quantum Set) 


(Classical Set) 


a = Kc{0-\qs) 


Mkc = Mk = 0.497 








a = kIi- \qs) 


Mki = Mk* = 0.896 


2.189 


2.462 


f3 = U^+\uds) 


Mj^ = 0.982 


1.648 


1.648 



TABLE 11. QMM parameters. The cluster color gyromagnetic ratios are fixed by the CMC constraint in terms of the cluster 
masses [Ge\^]. We fix the color octet cluster masses to equal the corresponding flavor equivalent K and K* masses. The only 
free parameter is the flavor singlet color octet mass (in bold type) adjusted to reproduce the lightest pentaquark candidate 
(Xo(1800)) observed in the SPHINX data [18-20]. 



Name 


x^0{r\L,s) 


QMM Mass [GoV] 
Quantum Parameters 


QMM Mass [GoV] 
Classical Parameters 


^o(r) 

^3(1 + ) 


XK*iJ2 |1'2) 


1.801 
1.967 
1.999 
2.064 


1.801 
1.967 
1.972 
2.049 



TABLE III. The comparison between calculated QMM exotic masses [GeV] using quantum consistent and classical param- 
eters for the gyromagnetic ratios. Note the different parameter sets only affect the K* gyromagnetic ratio, so the two lightest 
exotics with (i^c-^c) substructure are unaffected. The S.H.O. oscillator frequency was fixed by the A'^*(1440) — A'^(938) mass 
splitting: lj = 0.251 GeV. 
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TABLE IV. Masses and widths which characterize the narrow resonances fits to the p(7, K)[A, E°] SAPHIR data [21]. Mx, 



are the masses, Fx^ the widths, and Cx^ the signal strengths. N/S: not seen in this channel. N/F: peak not fitted. C.L.: 
confidence level. See text for details. 



JLab Data Fit 


p{e,e'K+)A 


p(e,e'A'+)E° 


Mxo[MeV] 


N/F 


O/A 


Fxo [MeV] 


N/F 


O/A 


Cxo 


N/F 


O/A 


C.L. 






Mat. [MeV] 


1896.5 ± 16.4 


N/F 


Fjv- [MeV] 


239.09 ± 121.4 


N/F 


Cat* 


1.18 ±0.01 


N/F 


C.L. 


> 99% 




Mxi [MeV] 


N/F 


1946.5 ± 4.2 


Fxi [MeV] 


N/F 


43.6 ± 16.0 


Cx, 


N/F 


1.53 ± 0.26 


C.L. 




> 99% 


Mx2 [MeV] 


N/F 


2001.7 ±5.2 


Tx2 [MeV] 


N/F 


54.1 ± 23.2 


Cx2 


N/F 


1.56 ±0.29 


C.L. 




> 99% 


[MeV] 


O/A 


O/A 


Fx3 [MeV] 


O/A 


O/A 


Cxs 


O/A 


O/A 


C.L. 







TABLE V. Masses and widths characteristics of the narrow resonances fit to the p{e, e'K) [A, E°] Jefferson Lab (JLab) 
data [25,26]. N/F: pealc not fitted. O/A: out of acceptance. C.L.: confidence level. Note the large width associated with the 
assumed N* peak (absent from the QMM spectrum). 
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SAPHIR 


JLab 












p{e, e K^)A 


p{e, e K^jX, 


/ T7r\ 




(1812.0 ±7.0) 


N/S 


N/S 


N/S 


N/S 


1801 


Vxq 


(56 ± 1) 


N/S 


N/S 


N/S 


N/S 


, 1 -I- > 


Cxo 


- 


N/S 


N/S 


N/S 


N/S 




Mn' 


N/S 


(1898.4 ±9.1) 


N/S 


(1896.5 ± 16.4) 


N/S 


Missing 


Vn' 


N/S 


(246.2 ± 57.0) 


N/S 


(239.1 ± 121.4) 


N/S 


(?■) 


Cn* 


" 


1.11 ± 0.04 


N/S 


1.18 ± 0.01 


N/S 




Mxi 


N/S 


N/S 


N/S 


(1946.5 ± 4.2) 


(1938.4 ± 6.7) 


1967 




N/S 


N/S 


N/S 


(43.6 ± 16.0) 


(239.9 ± 57.9) 


(1 ) 


Cxi 


" 


N/S 


N/S 


1.53 ± 0.26 


1.43 ± 0.03 




Mx2 


(1997.0 ± 7.0) 


N/S 


(1974.2 ± 3.2) 


O/A 


(2001.7 ±5.2) 


1999 




(91 ± 17) 


N/S 


(114.5 ±11.3) 


0/A 


(54.1 ± 23.2) 




CX2 




N/S 


1.10 ±0.04 


O/A 


1.56 ±0.29 




Mx:, 


(2052.0 ± 6) 


(2068.6 ± 4.3) 


(2072.7 ± 4.7) 


O/A 


O/A 


2064 


Tx3 


(35.0 ± 29) 


(112.1 ± 24.3) 


(113.5 ± 30.9) 


O/A 


O/A 


(1+) 






1.36 ±0.11 


1.27 ± 0.08 


O/A 


O/A 





TABLE VI. Experimental masses and widths of the narrow resonances fit to the ^^C{p,K~'')X SPHINX data [18-20], 
p(7,iC)[A, E°] SAPHIR data [21], and p(e, e'if )[A, S°] Jefferson Lab (JLab) data [25,26]. Mx are the masses [MeV], Fx the 
FWHMs [McV], and Cx the Breit-Wigner coupling strengths. N/S: not seen in this channel, O/A: out of acceptance. The 
QMM masses were calculated using the quantum consistent parameter set. 
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FIG. 1. The center-of-mass orbital angular momentum of the pentaquark system. 



a 











a p 

FIG. 2. Feynman diagram of the one gluon exchange between the the color octet clusters. 
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FIG. 3. A (left) and E'' (right) invariant mass distributions in photo-production from [21] compared to the model prediction 
of [2,3] when the narrow baryon X-resonances are taken into account (solid line), and without X (dot-dashed line). The dotted 
line in the A channel is the prediction from [27]. See text for details. 
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FIG. 4. Missing mass distribution for the kaon electro-production reaction [25,26]. Random coincidences and target wall 
contributions are subtracted. 
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FIG. 5. Invariant mass distributions of the A (loft) and (right) hypcrons. The model prediction is compared to the 
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the prediction of [2,3]. See text for details. 
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